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PREFACE 


This  report  was  prepared  by  Professor  J.  K.  Mitchell  and  Messrs.  J.  M. 
Chatoian  and  G.  C.  Carpenter  under  Contract  DACA  39-75-M0260  as  art  of  on- 
going work  conducted  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  under  the  sponsorship  of  the  Directorate  of  Military  Construction, 
Office,  Chief  of  Engineers,  under  Research  in  Military  Construction,  Project 
l»Al6ll02AT22,  Task  A2,  Work  Unit  005,  Study  of  Engineering  Classification 
of  Cohesionless  Soils. 

The  work  was  performed  under  the  direction  of  Drs.  J.  P.  Mulilis  and 
F.  C.  Townsend,  Research  Civil  Engineers,  Soil  Mechanics  Division,  Soils 
and  Pavements  Laboratory  (S&PLJ.  General  guidance  was  provided  by  the 
following  S&PL  personnel:  Messrs.  J.  P.  Sale  and  S.  J.  Johnson,  Chief  and 

Special  Assistant,  respectively,  S8PL;  and  Mr.  C.  L.  McAnear,  Chief,  Soil 
Mechanics  Division. 

Director  of  WES  during  this  study  and  preparation  and  publication 
of  this  report  was  COL  G.  H.  Hilt,  CE.  Technical  Director  was  Mr.  F.  R. 
Brown. 
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I . INTRODUCTION 


In  recent  studies  (Ladd,  1974;  Mulilis  et  al,  1975)  it  was  found 
that  the  method  of  sample  preparation  had  a significant  influence  on  the 
liquefaction  behavior  of  sands  under  cyclic  loading.  Figure  1 shows  the 
relationship  between  cyclic  deviator  stress  in  undrained  triaxial  tests 
as  a function  of  number  of  cycles  to  cause  initial  liquefaction.  Initial 
liquefaction  is  defined  as  the  condition  when  the  pore  pressure  has 
increased  to  a value  equal  to  the  confining  pressure. 

Differences  in  strength  under  cyclic  loading  of  the  magnitude  shown 
in  Fig.  1 have  important  implications  in  engineering  practice  if  recon- 
stituted samples  of  sand  are  to  be  used  for  evaluation  of  strengths  of 
sand  deposits  or  sand  fills  in  the  field.  In  view  of  the  difficulty  of 
obtaining  undisturbed  samples  of  sand  from  the  field  and  the  uncertainty 
as  to  how  undisturbed  "undisturbed  samples  ’ really  are,  it  is  nonetheless 
desirable  to  use  reconstituted  samples  if  possible. 

One  approach  to  the  selection  of  an  appropriate  laboratory  sample 
preparation  procedure  is  to  study  the  fabric  of  the  particulate  sand 
structure  formed  by  different  methods  in  conjunction  with  the  observed 
differences  in  liquefaction  behavior.  Then,  if  the  fabrics  of  field 
samples  can  be  evaluated,  the  reconstitution  of  samples  with  appropriate 
fabric  in  the  laboratory  should  be  possible.  With  the  exception  of  the 
results  of  studies  by  sedimentplogists  of  long  axis  orientations,  fabric 
details  of  sands  in  the  field  remain  largely  undetermined  at  the  present 
time.  Their  determination  should  be  relatively  easy.  In  addition  studies 
of  sand  fabric  are  useful  to  provide  a better  understanding  of  the  basic 
factors  controlling  the  engineering  properties  of  soils. 

Previous  work  by  Oda  (1972),  Mahmood  (1973),  and  Mahmood  and  Mitchell 
(1974)  has  established  that  the  static  strength  and  compressibility  behavior 


of  sands  and  silts  depends  on  the  fabric.  In  an  investigation  reported 
by  Mulilis  et  al  (1975)  it  was  shown  that  the  fabrics  of  samples  of 
Monterey  No.  0 sand  prepared  by  two  different  methods,  dry  pluviatior. 
and  moist  tamping,  were  different,  as  were  also  the  undrained  strengths 
in  undrained  cyclic  triaxial  tests.  Of  interest  in  their  study  were  the 
orientations  of  particle  long  axes  and  interparticle  contact  planes.  Their 
results  were  limited,  however,  in  that  insufficient  data  were  obtained 
to  obtain  a complete  distribution  of  contact  orientations,  and  only  one 
relative  density  was  studied.  Thus  the  conclusions  concerning  fabric 
effects  were  necessarily  tentative. 

The  purposes  of  the  present  study  were  to  make  a more  in-depth 
study  of  the  fabric  of  sand  samples  in  relation  to  method  of  compaction 
and  to  relate  these  fabrics  to  strengths  under  cyclic  loading.  In  addition 
behavior  in  drained  triaxial  compression  was  also  studied,  because  it  was 
considered  that  the  stress-strain  and  volume  change-strain  characteris- 
tics should  relate  to  both  the  fabric  and  the  liquefaction  behavior. 

All  testing  was  done  using  Monterey  No.  0 sand,  with  samples  pre- 
pared to  relative  densities  of  50  percent  and  80  percent.  Only  the 
behavior  at  a relative  density  of  50  percent  was  studied  by  Mulilis  et  al. 
Emphasis  in  the  fabric  study  was  on  grain  and  interparticle  contact 
orientations,  as  previous  work  by  Oda  (1972a,  1972b)  had  indicated  them 
to  be  of  particular  significance. 


II . SAMPLE  PREPARATION 


Introduction 

In  the  previous  study  by  Mulilis  et  al  (1975)  it  was  found  that  for 
Monterey  No.  0 sar.d  at  a relative  density  of  50  percent  the  greatest 
differences  in  strength  under  cyclic  loading  were  between  samples  pre- 
pared by  dry  pluviation  and  moist  vibration,  as  shown  in  Fig.  1.  Thus 
these  methods  were  adopted  for  this  investigation.  In  addition  it  was 
desired  to  compare  cyclic  strengths  and  fabrics  of  samples  of  the 
same  sand  at  a relative  density  of  80  percent.  As  samples  could  not  be 
prepared  to  this  higher  relative  density  using  moist  vibration,  a moist 
romping  procedure  was  used.  Samples  were  also  prepared  at  a relative 
density  of  50  percent  by  moist  tamping  for  comparative  purposes. 
Cylindrical  samples  prepared  by  all  methods  were  approximately  2.8  inches 
in  diameter  and  7.0  inches  high.  A split  mold  lined  with  an  impervious 
rubber  membrane  was  used. 

The  testing  program,  therefore,  encompassed  relative  densities  and 
compaction  methods  as  follows: 

Compaction  Method 
Dry  Pluviation 
Moist  Vibration 
Moist  Tamping 

Sand  Properties 

All  tests  were  done  using  Monterey  No.  0 sand,  a uniform  medium 
sand  with  rounded  to  subrounded  grains  consisting  predominantly  of  quartz 
with  some  feldspar  and  mica.  The  maximum  and  minimum  dry  densities  of 


Relative  Density 
50%  80% 

X X 

x 1 
X X 


this  sand  are  105.7  pcf  and  89.3  pcf,  respectively.  A grain  size  distri- 
bution curve  is  shown  in  Fig.  2 and  a grain  shape  histogram  of  particle 
length  to  width  ratios  is  shown  in  Fig.  3. 

Pluvial  Compaction 

Dry  sand  was  poured  from  a 1000  ml  flask  through  a one-hole  rubber 
stopper  containing  a nozzle  of  either  0.27  inch,  0.20  inch,  or  0.15  inch 
diameter.  Sand  free  fall  to  the  surface  of  the  sample  was  about  7 inches. 
The  relative  density  obtained  was  a function  of  the  intensity  of  particle 
rain,  which,  in  turn,  depended  on  nozzle  size  and  rate  of  flask  rotation 
above  the  sample  as  follows: 

a.  Decreased  nozzle  diameter  resulted  in  increased  relative  density. 

b.  Increased  flask  rotation  rate  resulted  in  increased  relative 
density. 

Compaction  by  Horizontal  High  Frequency  Vibrations 

The  sand  was  uniformly  mixed  to  an  initial  moisture  content  of 
8 percent.  Samples  were  fabricated  in  seven  one-inch  thick  layers,  each 
vibrated  under  a uniform  surcharge  pressure  of  0.8  psi.  Circumferential 
vibrations  were  applied  to  the  outside  of  the  mold  using  a BVI  Vibro-Graver 
at  a frequency  of  120  cps. 

To  prevent  overdens ificat ion  of  the  lower  layers  during  vibratioh  of 
the  upper  layers  the  "undercompaction"  procedure  described  by  Mulilis  et  al 
(1975)  was  used.  The  initial  as-compacted  relative  densities  of  each  one- 
inch  layer  to  give  a final  uniform  relative  density  of  50  percent  were: 


1* 


2, 


Relative  Density 
(%) 

47 

48 

49 

50 

51 

52 

53 

The  surface  of  each  layer  was  scarified  prior  to  placement  of  the 
next  layer.  The  sample  mixing  container  and  mold  were  filled  with  C02 
initially  to  facilitate  subsequent  saturation  for  drained  triaxial  and 
cyclic  load  testing. 

Compaction  by  Moist  Tamping 

The  hand  tamper  and  procedure  described  by  Mulilis  et  al  (1975), 
Fig.  4-2  and  p.  57,  were  used.  The  samples  were  prepared  in  seven  one- 
inch  thick  layers.  For  samples  at  50  percent  relative  density  the  same 
undercompaction  was  used  as  for  samples  prepared  by  mcist  vibration.  To 
obtain  uniform  samples  at  a relative  density  of  80  percent  the  following 
initial  layer  relative  densities  were  used: 


Layer  No. 
(from  bottom  up) 

1 

2 

3 

4 

5 

6 
7 
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Layer  No. 
(from  bottom  up) 

1 

2 

3 

4 

5 

6 
7 


Relative  Density 

(») 

77 

78 

79 

80 
81 
82 
83 
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The  top  surface  of  each  layer  was  scarified  prior  to  placement  of  the  next 
layer. 

Resin  Impregnation 

Moist  samples  were  first  dried  in  c forced  air  oven  at  60°C.  Warm 
air  was  drawn  slowly  through  the  sample  by  connecting  the  bottom  of  the 
sample  mold  to  a vacuum  of  about  1 psi  and  allowing  warm  air  to  enter  through 
the  top.  From  two  to  five  hours  were  required  to  complete  the  drying 
process. 

The  sample  to  be  impregnated  was  then  cooled  and  subjected  to  a 
vacuum  of  5 to  7 psi.  Koppers  polyester  resin  No.  1086-4,  mixed  with 
20  percent  (by  volume)  styrene  (Koppers  No.  VM-901)  to  reduce  viscosity 
and  1 percent  (by  volume)  organic  peroxide  (Witco  Chemical  HI-Point-180) 
as  catalyst,  was  introduced  at  the  bottom  of  the  sample.  After  saturation 
the  sample  was  left  undisturbed  until  initial  set  of  the  resin  had  occurred. 
The  sample  was  then  removed  from  the  mold  and  placed  in  a forced  air  oven 
at  60°C  for  24  hours  to  permit  the  resin  to  cure. 

Four  impregnated  samples  were  prepared  for  each  relative  density  and 
compaction  method,  for  a total  of  20  samples. 

Thin  Sections 

Three  thin  sections  were  cu*:  from  the  central  portions  of  three  of 
the  four  impregnated  samples  at  orientations  as  shown  in  Fig.  4.  By 
determining  grain  and  contact  orientations  in  each  of  these  thin  sections, 
the  distributions  in  all  directions  within  the  sample  could  be  determined. 
This  was  not  possible  in  the  preliminary  fabric  study  described  in  Mulilis 
et  al  (1975). 
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Petrographic  thin  sections,  about  60  pm  thick,  20  mm  wide,  and 
30  mm  long,  were  prepared  using  standard  methods. 

Photomicrographs  of  typical  thin  sections  of  samples  prepared  by 
each  of  the  methods  and  relative  densities  used  in  this  study  are  shown 
in  Figs.  5 through  9.  Vertical  thin  sections  (sections  parallel  to 
cylindrical  sample  axis)  are  shown  in  each  case. 

X-Radiographs 

Longitudinal  slices  about  3 mm  thick  were  cut  from  the  remaining 
cylindrical  sample  of  each  set  of  four.  Transmission  X-radiographs  were 
obtained  using  a medical  X-ray  machine.  X-radiographs  for  samples  pre- 
pared to  a relative  density  of  50  percent  are  shown  in  Fig.  10,  and  to  a 
relative  density  of  80  percent  in  Fig.  11. 


III.  FABRIC  ANALYSIS 


Fabric  elements  studied  quantitatively  were  the  three-dimensional 
distributions  of  normals  to  interparticle  contact  planes  and  the 
orientations  of  apparent  long  axes  in  the  plane  of  the  thin  section. 

Other  fabric  elements  that  may  have  ar,  important  influence  on  the  cyclic 
load  strength  and  other  mechanica1  properties,  such  as  pore  sizes  and 
pore  size  distributions,  particle  clustering,  if  any,  void  variability 
within  samples,  and  particle  size  segregation,  could  not  be  investigated 
within  the  limits  of  time  and  resources  available  for  this  study. 

Universal  Stage  Microscope 

A microscope  equipped  with  a universal  stage  was  used  to  determine 
particle  and  contact  orientations.  The  universal  stage,  shown  schematical- 
ly in  Fig.  12,  allows  independent  rotation  of  the  thin  section  about  three 
orthogonal  reference  axes.  The  thin  section  is  covered  by  the  upper  glass 
hemisphere,  H in  Fig.  12,  set  in  the  metal  mount  M.  The  thin  section  is 
aligned  and  moved  with  the  aid  of  the  sledge  R which  moves  vertically  in  a 
slot  cut  in  the  hemisphere  mount.  Rotation  can  be  made  about  the  NS,  EW, 
and  vertical  axes,  and  the  amount  of  rotation  is  measured  on  toe  appropriate 
vernier. 

The  use  of  the  universal  si  a tie  microscope  is  desi  ribed  in  detail  by 
Turner  and  Weiss  (1963)  . 

Sampling  Methods 

The  long  axis  orientations  and  contact  plane  orientations  can  be 
considered  as  two-dimensional  sample  populations  because  the  thickness  of 
the  thin  sections  is  very  small  compared  to  their  length  and  width.  Several 
methods  of  sampling  two-dimensional  populations  are  available  (Koch  and  Link, 
1970}  Miller  and  Kahn,  1962;  Kellerhals  and  Bray,  1971;  Pincus,  1953): 
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1.  Sample  the  entire  population. 

This  is  the  only  method  that  can  describe  the  true  orientation 
of  the  entire  population.  This  method  is  impractical  for  very  large 
populations,  as  is  the  case  here. 

2.  Grid  sampling. 

A grid  is  established  over  the  sample  and  only  those  grains  or 
contacts  immediately  below  the  grid  points  are  counted.  A variation 
of  this  method  is  to  assign  a number  to  each  grid  block  and  to  count 
entire  populations  within  randomly  selected  grid  blocks. 

3.  Transect  sampling. 

Every  member  of  the  parent  population  that  falls  under  a straight 
line  or  narrow  strip  across  the  sample  is  counted. 

4.  Random  sampling. 

Randomly  chosen  members  of  the  parent  population  are  measured. 

This  method  is  the  most  susceptible  to  sampling  bias. 

The  transect  sampling  method  was  used  in  this  investigation,  because 
(1)  the  population  was  too  large  to  sample  in  its  entirety,  (2)  no  grid 
attachment  was  available  for  the  microscope  and  (3)  it  was  necessary  to  avoid 
sampling  bias.  Traverses  about  eight  or  nine  grains  wide  and  about  one  inch 
long  were  chosen  parallel  to  the  long  axis  of  the  thin  sections.  Every  grain 
or  contact  within  this  zone  was  counted. 

Apparent  Longest  Axis  Orientations 

The  orientation  of  a particle  is  represented  by  the  inclination  of 
its  true  long  or  short  axis  with  respect  to  fixed  reference  axes.  Deter- 
mining the  true  orientation  of  the  long  or  short  axis  of  small  particles 
is  difficult.  In  practice  it  is  usually  done  by  analyzing  the  orientation 
of  apparent  longest  axes  as  seen  in  thin-seotions.  The  three-dimensional 
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orientation  of  grains  can  be  estimated  from  the  orientations  observed  in 
sections  cut  in  different  directions. 

In  this  study  the  orientations  of  apparent  longest  axes  were  deter- 
mined ir.  two  vertical  and  one  horizontal  section  for  each  sample. 
Orientations  were  determined  for  about  220  grains  in  each  thin  section 
by  measuring  the  angle  0^  between  the  apparent  longest  axes  and  reference 
axis  parallel  to  the  long  axis  of  the  thin-section  (Figure  4) . Both 
percent  frequency  histograms  (Figures  20  and  21)  and  rose  diagrams  (Figures 
22  and  23)  were  constructed  from  the  data  thus  obtained  using  an  azimuth 
class  interval  of  0 = 10°  to  show  the  concentrations  of  long  axes  in 
different  directions. 

Statistical  analyses  of  the  data  were  done  using  a vector  method  and 
test  of  significance  presented  by  Curry  (1956) . Uith  this  method  each 
long  axis  observation  is  considered  a vector  in  the  measured  direction  of 
unit  magnitude.  No  distinction  is  made  of  one  end  of  the  particle  from 
the  other.  All  measurements  are  therefore  made  in  the  range  0°  ^ 0 £ 180°, 
with  0=0°,  180°  corresponding  to  particle  axes  parallel  to  the  longi- 
tudinal axis  of  the  cylindrical  sample. 

A resultant  vector  calculated  on  the  180°  distribution  may  not 
reflect  the  true  central  tendency  of  the  distribution.  To  illustrate 
this  point  consider  the  following  example:  if  the  180°  range  lies  in  the 

eastern  semicircle,  the  distribution  has  no  west  components  at  all.  North 
components  would  tend  to  cancel  south  components  so  the  resultant  vector 
would  always  possess  a strong  easterly  central  tendency  even  if  the  true 
central  tendency  lay  close  to  the  north-south  line.  To  remedy  this  the 
angles  of  the  observation  vectors  are  doubled  before  computing  the 
components,  thus  obtaining  a nonsymmetric  periodic  distribution. 
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The  calculations  are: 


N-S  component  = In  cos20 
E-W  component  = In  sin20 


where  n = number  of  observations  oriented  at  0. 


_ £n  sin20 
tan  20  - In  cosTG" 


vector  magnitude. 


vector  direction, 


r = /(In  sin20)  2 + (In  cos20): 


L = ■=—  x 100% 
In 


1 .In  sin20 
0 = 2 arctan  In"  "cos 20 


The  vector  direction  is  a measure  of  the  preferred  orientation 
direction  of  long  axes  of  sand  grains.  It  has  the  advantage  of  being 
independent  of  the  choice  of  origin  (Curry,  1956). 

The  vector  magnitude  varies  from  0 percent  to  100  percent.  A random 
(i.e.  : uniform)  distribution  of  orientations  will  give  a vector  magnitude 
of  0 percent.  A vector  magnitude  of  100  percent  means  that  all  orientations 
are  exactly  the  same  or.  in  the  case  of  grouped  data,  all  orientations  lie 
within  the  same  class  interval.  According  to  Curry  (1956)  L is  a sensitive 
measure  of  dispersion  and  is  comparable  to  standard  deviation  but  has  the 
advantage  of  being  independent  of  the  choice  of  origin. 

Curry  (1956)  also  adapted  the  vector  magnitude,  L,  for  use  in  the 
Rayleigh  test  of  significance.  Rayleigh's  equation  is: 

-r2/n 

P = e (5) 

where  p = probability  of  obtaining  a greater  amplitude  by  pure  chance 
combinations  of  random  phases, 
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r - resultant  an^litude  obtained, 
n = number  of  observations, 
e - base  of  natural  logarithms. 

For  the  resultant  amplitude,  r expressed  in  toms  of  the  vector  magnitude. 


then 


r 


Ln 

100 


(-L2n  x 10“4) 

p = e 


(6) 

(7) 


Any  level  of  significance,  p can  be  used.  In  practice  the  0.05  level 
is  commonly  used  for  the  study  of  sand-grain  orientations.  The  0.05  level 
means  that  there  are  only  5 chances  in  100  of  the  determined  distribution 
being  due  to  chance. 

It  was  assumed  that  to  establish  randomness  there  Jhould  be  a minimum 
of  10  observations  per  10°  azimuth  class  interval.  There  were  18  azimuthal 
intervals  thereby  gjving  a minimum  sample  of  180  observations.  After 
180-200  grains  were  sampled  for  a thin  section  the  data  were  tabulated  and 
L,  0,  and  p determined.  If  the  significance  level  was  less  than  0.05 
(i.e.:  p > 0.05)  additional  observations  were  made  until  p < 0.05. 

One  test  was  conducted  to  determine  the  effect  of  sample  size  using 
the  same  thin  section.  Two  groups  of  data  were  obtained,  one  consisting 
of  198  observations  and  the  other  of  484  observations.  The  determined  dis- 
tributions were  very  similar,  and  the  mean  vector  directions,  0,  were  only 
10°  apart.  This  is  considered  to  be  close  agreement  as  any  one  observation 
is  only  accurate  to  ±2°  to  4°,  depending  upon  the  grain  shape.  Both  tests 
gave  Rayleigh  levels  of  significances  greater  than  0.05.  It  was  therefore 
concluded  that  a sample  consisting  of  about  200  observations  would  provide 
as  accurate  a measure  of  long  axis  orientation  distribution  as  one  consist- 
ing of  over  450  observations,  so  long  as  p £ 0.05. 
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Intarparticle  Contact  Orientations 


The  orientation  of  any  interparticle  contact  can  be  represented  by 
the  normals,  1 , N^2,  to  the  tangent  plane,  as  shown  in  Fig.  13.  A method 
for  analyzing  the  orientations  of  contact  plane  normals  has  been  developed 
by  Oda  (1972)  and  was  used  also  for  this  investigation. 

It  is  convenient  to  express  vector  orientations  in  terms  of  spherical 
coordinates.  Spherical  and  Cartesian  coordinates  are  related  according  to 
x = r sin  2 cos  a (8a) 

y = r sin  2 sin  a (8b) 

z = r cos  2 (8c) 


where  r is  radius  or  vector  length,  and  angles  a and  2 describe  vector 
direction  (Fig.  13) . 

An  element  dA  of  the  surface  area  of  a sphere  is  given  by 
dA  Q = rd2’r  sin2  da  = r2  sin2  dad2  (9) 


for  a sphere  of  unit  radius 

dA  0 = sin2  dad 3 
a,p 


Integration  yields 


sin2  da  d2 


(10) 


(11) 


Consider  M contacts  within  an  assemblage  of  granular  particles.  At 
each  interparticle  contact  there  are  two  contact  surfaces,  one  belonging 
to  grain  1 and  the  other  to  grain  2.  Normal  directions  N^1  and  N^2,  Figure 
13,  are  perpendicular  to  the  tangent  plane  at  the  contact  point.  The 
direction  of  N^2  is  described  by  angles  a and  2;  whereas,  the  direction  of 
NL  1 is  described  by  a+m  and  tt-B. 
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The  distribution  of  interparticle  contact  normals  can  be  shown  by 
plotting  their  piercing  points  on  the  sphere  or  its  graphical  equivalent, 
the  equal  area  stereonet  (Turner  and  Weiss,  1963) . The  total  number  of 
plotted  points  is  2M,  twice  the  number  of  contacts. 

The  total  number  of  contact  normals  within  the  angular  ranges 
a+da  and  3+d8  is 

2M  E(a,B)  sinB  dadB  (12) 

where  E(a,3)  is  the  probability  density  of  points  of  contact  within  the 
angular  intervals,  a function  of  both  o and  8.  All  plotted  points  will  be 
included  if  (12)  is  integrated  over  a = 0°  to  360°  and  8 = 0°  to  180'°. 

Thus  E(a,8)  must  satisfy 

E(ot,B)  sin8  dadB  = 1 (13) 


In  the  present  investigation  E(0C,B)  is  independent  of  a because  of 
the  axisymmetric  character  of  the  samples  studied.  Thus  equation  (13) 
becomes 

2 it  f E (8)  sinB  dB  = 1 (14) 

•/3 


As  vectors  N.  1 and  N.2  are  colinear  E ( 8 ) equals  E(tt-8).  For  example, 
ri 

the  probability  density  in  the  region  8 = 0°  to  8 = 10°  is  the  same  as  in 
the  region  8 = 180°  to  8 = 170°.  As  a result  only  the  vector  orientations 
in  one  hemisphere  need  be  measured  and  plotted  and  the  results  doubled  to 
obtain  the  density  distribution  over  the  entire  sphere.  Thus  equation  (14) 
becomes 


2 [2tt 


E(8)  sinB  dB] 


1 


(15) 


M = 41TM 


E(8) 


sinBdB  + + 


E(3)  sinB  dB 


(16) 


lU 


or 


then 


If  is  the  total  number  of  normals  in  the  range  8 = 0°  to  0 = 10°, 


M 


1 


E (0)  sin$d8 


(17) 


Integration  gives  the  mean  value  of  E^fB)  within  the  range  0°  £ 8 <_  10°, 

M1 

E (0)  = ± (18) 

4isM(cos  0°  cos  10°) 

Mean  values  of  E(0)  for  other  intervals  are  calculated  in  a similar 
manner. 

For  an  isotropic  fabric  E^(8)  = constant  = . Values  of  E^(0)> 

mean  a greater  than  average  proportion  of  contact  normals  at  inclination 
8^;  wheres  E^(8)  < ^ indicates  a less  than  average  proportion  of  con- 

tact normals  at  inclination  8.. 

l 
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IV.  LIQUEFACTION  BEHAVIOR 


Cyclic  triaxial  tests  were  done  following  the  procedures  used  by 
Mulilis  et  al  (1975).  The  purpose  of  these  tests  on  samples  prepared  to 
a relative  density  of  50  percent  was  to  confirm  that  our  sample  prepara- 
tion techniques  gave  results  comparable  to  those  of  Mulilis  et  al.  If  so, 
then  it  seemed  reasonable  that  the  fabric  measurements  and  drained  tri- 
axial test  results  could  be  applied  to  the  findings  of  the  earlier  study 
as  well.  Several  tests  were  done  by  each  method,  and  the  results  agreed 
with  the  curves  shown  in  Fig.  14. 

Additional  cyclic  load  triaxial  tests  were  done  on  samples  at  a 
relative  density  of  80  percent.  The  results  are  shown  in  Fig.  15.  It 
may  be  seen  that  even  at  this  high  relative  density  samples  prepared  by 
moist  tamping  are  much  stronger  than  samples  prepared  by  dry  pluviation. 

It  is  also  apparent  that  nozzle  size  for  preparation  of  samples  by  dry 
pluviation  has  a consistent  and  significant  effect,  with  the  use  of  the 
smaller  diameter  nozzle  (0.20  inch)  giving  stronger  samples. 
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V.  DRAINED  TRIAXIAL  COMPRESSION  BEHAVIOR 


Constant  rate  of  strain  drained  triaxial  compression  tests  were 
done  on  saturated  samples  prepared  at  50  percent  and  80  percent  relative 
density  by  the  different  methods.  The  purpose  of  these  tests  was  to 
determine  if  consistent  relationships  existed  between  the  static  stress- 
strain  and  volume  change  behavior,  the  fabric,  and  the  liquefaction 
behavior. 

After  compaction  the  air  in  the  voids  of  the  samples  was  replaced 
by  C02  to  facilitate  saturation.  A back  pressure  of  about  30  psi  was 
required  to  attain  essentially  complete  saturation  (B  > 0.94).  All 
samples  were  equilibrated  under  an  effective  isotropic  confining  pressure 
0^  of  8 psi.  Drained  triaxial  compression  was  carried  out  using  an 
arbitrarily  selected  constant  rate  of  axial  deformation  of  0.16  mm/min. 
This  corresponded  to  about  0.09  percent  axial  strain  per  minute.  Volume 
changes  were  measured  to  the  nearest  0.01  cc. 

Behavior  at  50  Percent  Relative  Density 

Deviator  stress  vs.  axial  strain  and  volumetric  strain  vs.  axial 
strain  for  each  of  two  tests  on  samples  at  50  percent  relative  density 
prepared  by  dry  pluviation,  moist  tamping  and  moist  vibration  are  shown 
in  Fig.  16.  Averaged  curves  for  the  two  tests  by  each  method  are  shown 
in  Fig.  17. 

It  may  be  seen  that  the  stiffness  at  low  strains,  as  indicated  by 
the  slopes  of  the  stress-strain  curves,  increases  in  the  order  pluviated  < 
moist  tamped  < moist  vibrated.  Also  the  amount  of  volumetric  compression 
preceding  dilation  increases  in  the  order  moist  vibrated  < moist  tamped  < 
pluviated.  Both  of  these  findings  are  consistent  with  the  differences  in 
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undrained  strength  in  cyclic  loading  shown  in  Fig.  14.  The  lower  the 
modulus  and  the  higher  the  tendency  for  volumetric  compression,  the 
greater  the  rate  of  stress  transfer  to  the  pore  water  under  undrained 
conditions,  and  therefore  the  fewer  cycles  of  loading  needed  to  cause 
failure. 

That  different  methods  of  sample  preparation  result  in  different 

sand  fabrics  and  differences  in  the  stress-strain  and  volume  change 

strain  behavior  has  been  well  established  in  previous  studies;  e.g.,  Oda 

(1972),  Arthur  and  Menzies  (1972,  Mahmood  (1973),  Mahmood  and  Mitchell  (1974). 

From  Fig.  14  it  may  be  seen  that  cyclic  stress  ratios  0^/20  ' required 

to  cause  liquefaction  after  10  to  100  cycles  range  from  about  0.3  to  0.5.  For 

the  value  of  O'  =8  psi  (55.2  kN/m2)  used  in  these  tests  the  corresponding  cyclic 
o 

deviator  stress  is  ±4.8  to  8 psi  (33  to  55  kN/m2).  Fig.  17  shows  that  the 
stiffnesses  and  volume  changes  associated  with  drained  uniaxial  compression 
under  stresses  of  these  magnitudes  are  essentially  the  same  for  all  three 
compaction  methods.  Thus  for  the  first  application  of  +occj  in  a cyclic  load 
test  there  should  be  essentially  no  difference  in  behavior  for  the  different 
samples. 

It  is  at  higher  values  of  deviator  stress  that  the  curves  begin  to 

diverge  in  the  drained  tests.  In  the  undrained  cyclic  load  tests  the 

divergence  must  begin  with  the  stress  reversal  (-0  ) . In  this  phase 

cd 

the  cell  pressure  becomes  0^  and  the  axial  stress  is  O^,  causing  a reversal 
in  shear  stress  direction.  Thus  the  uniaxial  compression  ana  cyclic 
load  tests  are  no  longer  strictly  comparable.  Nonetheless,  a distinct 
correlation  between  the  stiffness  and  volume  change  characteristics  in 
drained  compression  and  resistance  to  cyclic  loading  is  evident  for  samples 
at  50  percent  relative  density. 
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Behavior  at  80  Percent  Relative  Densit; 

Stress  vs.  strain  and  volume  change  vs.  strain  data  for  samples 
prepared  to  a relative  density  of  80  percent  by  dry  pluviation  and  moist 
tamping  are  shown  in  Fig.  18.  Although  the  peak  strength  and  the  dilation 
are  greater  for  the  moist  tamped  samples,  the  tangent  moduli  for  stresses 
up  to  about  50  percent  of  the  strength  and  the  volumetric  compressions 
prior  to  dilation  are  essentially  the  same  for  samples  prepared  by  both 
methods . 

By  analogy  with  the  results  at  50  percent  relative  density  (Figs.  14, 
16,  and  17)  it  might  be  anticipated  that  differences  in  cyclic  load 
strength  for  the  samples  at  80  percent  relative  density  would  be  relatively 
small.  Fig.  15  shows  that  this  is  not  the  case. 

Differences  in  drained  compression  behavior  were  observed 
for  samples  prepared  by  pluvial  compaction  using  different  nozzle  sizes, 
as  may  be  seen  in  Fig.  19.  Although  the  undrained  cyclic  loading  strength 
was  also  influenced  by  nozzle  size,  Fig.  15,  the  differences  were  nowhere 
near  as  great  as  between  the  pluviation  and  moist  tamping  compaction 
methods . 

Thus  at  high  relative  density  the  strength  under  cyclic  loading  does 
not  appear  to  be  as  distinctly  related  to  the  drained  triaxial  compression 
behavior  as  it  does  at  low  relative  density.  This  is  illustrated  further 
by  the  data  in  Table  1 which  compare  values  of  maximum  compressive 
volumetric  strain  preceding  dilation  and  secant  modulus  at  1 percent  axial 
strain. 
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Table  1 


Effect  of  Sample  Preparation  Method  on  Volume 
and  Stiffness  of  Monterey  No.  0 Sand 
in  Undrained  Triaxial  Compression 


Test 

No. 

Relative 
Density  (%) 

Preparation 

Method 

Maximum  Compressive 
Volumetric  Strain  (%) 

3,  4 

50.0 

Pluviation1 

0.105 

6,  8 

50.3 

Moist  Tamping 

0.051 

9,  10 

50.0 

Moist  Vibration 

0.044 

11,  12 

79.8 

Moist  Tamping 

0.050 

15,  16 

80.3 

Pluviation1 

0.049 

5,  7 

79.7 

2 

Pluviation 

0.036 

‘o^"  Dia.  Nozzle 
20.15"  Dia.  Nozzle 
Preceding  Dilation 

Values  are  the  average  of  two  tests 
= 8 psi  (55.2  kN/m2)  for  all  tests 


Change 

Secant  Modulus  at 
1%  Axial  Strain 
(kN/m2) 

140 

192 

210 

246 

207 

232 
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VI.  APPARENT  LONG  AXIS  ORIENTATIONS 


The  distributions  of  apparent  long  axis  orientations  as  seen  in  the 
thin  sections  and  the  directions  and  intensities  {vector  length,  L)  of 
preferred  orientations  were  determined  using  the  procedures  previously 
described  The  longitudinal  axis  (vertical  axis)  of  the  original  cylindri- 
cal sample  was  used  as  the  reference  axis,  for  tne  vertical  sections  and 
long  axis  orientations  are  expressed  in  terms  of  0,  the  angle  between  the 
reference  axis  and  the  long  axis.  Thus  a value  of  0 = 0°,  180°  means 
particles  are  oriented  in  a vertical  direction;  whereas,  0 = 90°,  270° 
indicates  horizontal  (radial)  orientation. 

The  choice  of  x and  y axes  for  referencing  of  long  axes  in  the 
horizontal  sections  was  arbitrary.  Because  of  radial  symmetry  preferred 
orientations  in  horizontal  thin  sections  should  be  small  and  reflect  only 
local  inhomogeneities  caused  by  uncontrollable  variations  during  sample 
preparation.  This  expectation  of  randomness  is  borne  out  oy  the  test 
results  in  Table  2 which  list  for  each  horizontal  section  the  magnitude  of 
preferred  orientation  vector  L and  its  direction  8.  It  may  be  seen  that 
in  most  cases  L is  small,  and  there  is  no  predominant  direction  for  0 
associated  with  any  of  the  sample  preparation  methods. 

Results  for  the  vertical  sections  are  listed  in  Table  3.  It  may  be 
seen  that  the  differences  among  the  values  for  the  six  thin  sections  studied 
for  each  relative  density  and  compaction  method  are  small.  Thus  any  one 
of  the  sample  preparation  methods  used  would  appear  to  produce  reasonably 
reproducible  samples.  Average  values  for  the  orientation  intensity  vector 
L and  the  orientation  direction  are  given  in  Table  4. 

Frequency  histograms  and  orientation  rose  diagrams,  both  based  on  10° 
azimuth  class  intervals,  were  prepared  from  the  data  for  each  thin  section. 
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Table  2 


Long  Axis  Orientations  of  Sand  Particles  in 
Horizontal  Thin  Sections 


Compaction 

Method 

Nominal 
Relative 
Density  (%) 

Sample 

Number 

Number  of 
Observations 
n 

Vertor 
Length 
L (%) 

Vector 

Direction 

9 (°) 

Pluviation 

50 

P1H 

234 

2.5 

54 

P2H 

216 

2.8 

83 

P3H 

234 

1.9 

133 

80 

P4H 

240 

14.1 

24 

P5H 

234 

3.0 

124 

P6H 

234 

7.0 

37 

Moist  Tamped 

50 

MT1H 

240 

12.1 

7 

MT2H 

234 

6.8 

135 

MT3H 

216 

8.9 

104 

80 

MT4H 

240 

3.8 

74 

MT5H 

216 

4.3 

36 

MT6H 

216 

4.4 

114 

Moist  Vibrated 

50 

MV1H 

240 

3.8 

23 

MV2H 

234 

3.6 

91 

MV3H 

240 

2.2 

121 

*Y-Y  direction  in  Fig.  4 represents  0=0°. 
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Table  3 


Long  Axis  Orientations  of  Sand  Particles  in 
Vertical  Thin  Sections 


Compaction 

Method 

Nominal 
Relative 
Density  (%) 

Sample 

Number 

Number  of 
Observations 
n 

Vector 
Length 
L (%) 

Vector 
Direction 
6 (%) 

Pluviation 

50 

P1V1 

286 

18.3 

86 

P1V2 

234 

13.5 

85 

P2V1 

234 

19.2 

87 

P2V2 

234 

21.5 

97 

P3V1 

240 

19.3 

86 

P3V2 

216 

16.7 

95 

80 

P4V1 

240 

28.5 

99 

P4V2 

234 

38.2 

77 

P5V1 

240 

31.9 

101 

P5V2 

240 

28.2 

76 

P6V1 

216 

31.2 

101 

P6V2 

234 

36.6 

80 

Moist  Tamped 

50 

MT1V1 

234 

18.5 

86 

MT1V2 

216 

18.7 

101 

MT2V1 

234 

14.9 

84 

MT2V2 

234 

17.3 

101 

MT3V1 

234 

19.5 

83 

MT3V2 

234 

18.5 

90 
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Table  3 (Cont.) 


Compaction 

Method 

Nominal 
Relative 
Density  (%) 

80 

Sample 

Number 

Number  of 
Observations 
n 

Vector 
Length 
L (%) 

Vector 

Direction’ 

F <%) 

Moist  Tamped 

MT4V1 

234 

12.5 

86 

MT4V2 

216 

15.8 

78 

MT5V1 

216 

16.3 

84 

MT5V2 

216 

13.0 

76 

MT6V1 

234 

12.2 

82 

MT6V2 

234 

12.4 

100 

Moist  Vibrated 

50 

MV1V1 

234 

10.2 

78 

MV1V2 

240 

13.2 

105 

MV2V1 

234 

11.2 

71 

MV2V2 

240 

11.1 

106 

MV3V1 

234 

11.4 

77 

MV3V2 

240 

13.4 

100 

*Z-Z  direction  in  Fig.  4 represents  0-0°. 
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Table  4 


Summary  of  Long  Axis  Orientation  Data 
Monterey  No.  0 Sand 


Two  Vertical  Thin  Sections  as  Shown  in  Fig.  4 from 
Three  Samples  Prepared  by  Each  Method.  Between  216 
and  286  Observations  per  Thin  Section. 


Compaction 

Method 

Relative 
Density  {%) 

Vector 
Length 
L (%) 

Vector 

Direction 

Pluviation 

50 

18.1 

89.3 

80 

32.4 

89.0 

Moist  Tamped 

50 

17.9 

90.8 

80 

13.7 

84.3 

Moist  Vibrated 

50 

11.8 

89.5 
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Only  typical  examples  are  shown  here,  corresponding  to  the  different 
relative  densities  and  compaction  methods.  Figs.  20  and  21  are  histograms 
which  compare  samples  prepared  by  different  methods  to  relative  densities 
of  50  percent  and  80  percent,  respectively.  The  horizontal  dashed  line 
in  each  figure  represents  a completely  random  orientation.  The  rose 
diagrams  in  Figures  22  and  23  are  similar  comparisons.  The  light  circles 
shown  in  these  figures  represent  a random  orientation. 

The  vector  magnitudes  for  samples  prepared  by  moist  tamping  at  80 
percent  relative  density  and  by  moist  vibration  at  50  percent  relative 
density,  Tables  3 and  4,  are  small  and  are  close  to  the  lower  level  of 
significance  according  to  the  Rayleigh  test  described  previously.  On  the 
other  hand  the  vector  magnitudes  and  sample  populations  for  the  other 
thin  sections  give  a high  level  of  significance.  This  is  evidence  that 
the  moist  tamped  samples  at  80  percent  relative  density  and  the  moist 
vibrated  samples  at  50  parcent  relative  density  havt  a nearly  random 
distribution  of  long  axis  orientations.  These  were  the  strongest  samples 
in  both  undrained  cyclic  triaxial  loading  and  drained  triaxial  compression. 

Definite  preferred  orientation  of  long  axes  in  a direction  normal 
to  the  longitudinal  axis  was  developed  by  pluviation.  Significant 
orientation  in  the  same  direction  was  also  obtained  by  moist  tamping  to  a 
relative  density  of  50  percent. 

Thus,  for  this  sand  and  for  the  three  methods  of  compaction  used,  it 
appears  that  for  a given  relative  density  the  greater  the  intensity  of 
preferred  long  axis  orientation  in  a direction  normal  to  the  direction 
of  the  applied  cyclic  deviator  stress,  the  less  the  resistance  to  lique- 
faction under  triaxial  loading  conditions. 
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VII.  INTERPARTICLE  CONTACT  ORIENTATIONS 


Interparticle  contact  orientations,  expressed  in  terms  of  the  normal 
to  the  interparticle  contact  plane  both  by  equal  area  stereonets  and  by 
the  density  function  E{0),  developed  in  a previous  section,  were  deter- 
mined for  each  fabric  and  relative  density  combination  studied.  As  was 
the  case  for  the  distribution  of  long  axis  orientations,  the  distributions 
of  interparticle  contact  normals  were  remarkably  similar  for  the  three 
samples  prepared  by  each  method.  This  provides  further  support  for  the 
concept  that  different  reproducible  fabrics  are  formed  by  the  different 
compaction  methods. 

The  equal  area  stereonets  are  shown  in  Figs.  24  through  28.  On 
these  diagrams  blank  areas  indicate  contact  normal  directions  within  which 
less  than  one  percent  of  the  normals  intersect  one  percent  of  the  surface 
area  of  a sphere  containing  all  the  possible  directions.  Light  stippled 
areas  enclose  1 to  2 percent  of  the  contact  normals  per  one  percent  of 
the  circumscribing  sphere.  The  heavy  stippled  areas  indicate  more  than 
2 percent  of  the  normals  per  one  percent  of  the  sphere  surface.  Thus  the 
darker  the  area  the  greater  the  concentration  of  interparticle  contact 
normals. 

The  distribution  functions  E(B)  for  samples  prepared  by  pluvial 
compaction  to  relative  densities  of  50  percent  and  80  percent  are  compared 
in  Fig.  29.  Similar  plots  for  samples  prepared  by  moist  tamping  are  shown 
in  Fig.  30.  The  results  for  a sample  prepared  by  moist  vibration  to  a 
relative  density  of  50  percent  are  shown  in  Fig.  31.  In  these  figures 
the  light  horizontal  dashed  line  indicates  the  E($)  distribution  for  a 
completely  random  distribution. 


E(8)  distributions  for  samples  prepared  to  a relative  density  of 
50  percent  by  three  methods  are  compared  in  Fig.  32.  A comparison  of 
two  methods  at  a relative  density  of  80  percent  is  given  in  Fig.  33. 

The  stereonets  for  samples  prepared  by  pluviation,  Figs.  24  and  27, 
and  the  corresponding  E(8)  distributions,  Fig.  29,  show  that  interparticle 
contacts  are  relatively  randomly  oriented  at  50  percent  relative  density. 

At  80  percent  relative  density,  however,  there  is  a predominance  of  inter- 
particle contact  normals  oriented  at  0°  to  30°  of  the  vertical.  This 
corresponds  tc  a significant  increase  in  the  proportion  of  contact  planes 
oriented  at  0 to  30°  to  the  horizontal. 

Both  moist  tamping,  Figs.  25  and  30,  and  moist  vibration,  Figs.  26 
and  31,  show  concentrations  of  interparticle  normals  oriented  in  the  range 
of  0 to  40°  of  the  vertical  for  samples  at  50  percent  relative  density, 
indicating  preferred  orientation  of  contact  planes  at  0 to  40°  with  the 
horizontal.  A sLidlar  relationship  is  seen  for  samples  prepared  by 
moist  tamping  at  a relative  density  of  80  percent,  Figs.  28  and  30. 

Comparison  of  the  E(8)  distributions  for  the  three  compaction  methods 
and  samples  at  50  percent  relative  density.  Fig.  32,  shows  that  the  pro- 
portion of  contact  planes  oriented  in  the  range  of  0°  to  30°  with  the 
horizontal  increases  in  the  order  pluviation,  moist  tamping,  moist  vibra- 
tion. This  is  the  same  order  of  strength  increase  under  cyclic  loading. 

It  agrees  also  with  the  order  of  stiffness  increase  and  dilation  in 
drained  triaxial  compression.  It  appears,  therefore,  that  the  greater 
the  proportion  of  interparticle  contacts  oriented  normal  to  the  direction 
of  the  major  principal  stress  in  triaxial  compression  the  greater  the 
resistance  to  deformation  and  liquefaction  of  Monterey  No.  0 sand  at  a 
relative  density  of  50  percent. 
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At  80  percent  relative  density  samples  prepared  by  pluviation  have 
a somewhat  greater  proportion  of  interparticle  contacts  oriented  close  to 
the  horizontal  (0  to  30°),  Fig.  33,  than  do  samples  prepared  by  moist 
tamping.  This  result  was  somewhat  surprising  in  view  of  the  very  much 
higher  strength  under  cyclic  loading  for  sample  prepared  by  moist 
tamping,  Fig.  15.  The  differences  in  strength  and  volume  change  behavior 
in  drained  triaxial  compression,  Fig.  18,  were  relatively  small  by 
comparison,  but  they  were  also  in  a contrary  direction  to  the  strength- 
contact  orientation  relationship  found  at  50  percent  relative  density  for 
the  pluviation  and  moist  tamping  methods. 
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VIII.  SAMPLE  UNIFORMITY 


The  uniformity  of  samples  on  a macroscale  is  revealed  by  the  X- 
radiographs  of  vertica  . slices  shown  in  Fig.  10  for  a relative  density 
of  50  percent  and  Fig.  11  for  a relative  density  of  80  percent.  Both 
moist  tamping  and  moist  vibration  produce  narrow  (~1  mm)  bands  of  high 
density  at  the  layer  boundaries.  The  difference  between  the  density  in 
this  band  and  the  remainder  of  each  1 inch  layer  at  50  percent  relative 
density  is  greater  than  at  80  percent  relative  density,  as  may  be  seen 
by  comparing  the  X-ridiographs  for  moist  tamped  samples  in  Figs.  10  and 
11.  The  banding  is  more  intense  in  moist  vibrated  samples  than  in  moist 
tamped  samples,  Fig.  10. 

Thin  (1-2  mm) , wavy  layers  of  alternating  high  and  low  density  are 
observed  in  samples  prepared  by  pluviation.  Differences  in  these  layers 
are  less  at  80  percent  relative  density  than  at  50  percent  relative 
density.  Macroscale  fabric  differences  at  80  percent  relative  density 
resulting  from  pluvial  and  tamping  compaction  methods  are  less  pronounced 
than  at  50  percent  relative  density,  even  though  the  cyclic  load  strength 
differences  are  large  at  80  percent  relative  density. 

In  a previous  study  (Mulilis  et  al,  1975)  it  was  shown  that  the 
inclusion  of  layers  of  differing  density  within  a sample  could  not  be 
the  complete  cause  of  measured  differences  in  the  liquefaction  behavior 
of  samples  prepared  to  the  same  average  relative  density  by  different 
methods.  Nonetheless,  differences  of  the  type  shown  in  Figs.  10  and  11 
can  be  expected  to  influence  the  mechanical  properties.  The  full  magnitude 
of  this  influence  can  only  be  established  after  further  study. 

It  seems  reasonable  to  conclude  also  that  moist  tamping  and  moist 
vibration  cannot  reasonably  be  expected  to  produce  macro- fabrics  that 
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correspond  to  water  or  wind  sedimented  sand  deposits.  These  methods  may 
however,  model  in  some  way  the  fabrics  of  sand  fills  densified  in  layers 
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IX.  DISCUSSION 


The  results  presented  in  the  preceding  sections  have  shown  the 
following  relationships  for  Monterey  No.  0 sand. 

1.  The  fabrics  produced  by  each  of  the  sample  preparation  methods 
studied  were  quite  reproducible. 

2.  Samples  prepared  to  both  50  percent  and  80  percent  relative 
density  by  dry  pluviation  are  considerably  weaker  under 
undrained  cyclic  loading  than  samples  prepared  by  moist 
vibration  or  by  moist  tamping. 

3.  Samples  prepared  by  dry  pluviation  are  less  stiff,  compress 
more  prior  to  dilation  (at  50  percent  relative  density) , and 
dilate  less  than  samples  prepared  by  moist  vibration  or  moist 
tamping  when  tested  in  drained  triaxial  compression. 

4.  Differences  in  the  drained  stress-strain  and  volume  change- 
strain  behavior  are  considerably  less  at  80  percent  relative 
density  for  different  methods  of  compaction  than  are  differences 
in  the  liquefaction  behavior. 

5.  Preferred  orientations  of  particle  long  axes  develop  in  the 
horizontal  direction  for  each  of  the  methods  of  compaction 
studied. 

6.  The  intensity  of  preferred  long  axis  orientation,  as  measured 
by  vector  length  L in  Table  4,  increases  in  the  order  moist 
vibration,  moist  tamping,  dry  pluviation,  with  the  particles 
very  nearly  randomly  oriented  in  the  samples  prepared  by  moist 
vibration. 

7.  At  a given  relative  density  the  more  random  the  particle 
orientation  (lower  L)  the  higher  the  strength  under  undrained 
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cyclic  triaxial  loading  and  the  greater  the  stiffness  and 
dilation  in  drained  triaxial  compression. 

8.  At  a relative  density  of  50  percent  the  proportion  of  inter- 
particle contact  plane  orientations  in  the  range  0 to  ±30°  of 
the  horizontal  increases  in  the  order  pluviation,  moist  tamping, 
and  moist  vibration.  The  proportion  of  contact  planes  oriented 
in  these  directions  is  greater  than  at  other  orientations  for 
all  methods  of  compaction. 

9.  At  80  percent  relative  density  a greater  proportion  of  inter- 
particle contacts  is  oriented  at  0 to  ±30°  with  the  horizontal 
for  dry  pluviation  than  for  moist  tamping,  even  though  the 
strength  in  cyclic  loading  is  considerably  greater  for  samples 
prepared  by  moist  tamping. 

10.  Substantial  variations  in  sample  density  in  a longitudinal 

direction  could  be  seen  in  X-radiographs , and  these  variations 
differed  for  sample  preparation  by  the  different  methods. 

These  results  taken  collectively  indicate  that  the  fabric  elements 
represented  by  particle  long  axis  orientation  and  interparticle  contact 
plane  orientation  are  strongly  related  to  both  the  undrained  strength  under 
cyclic  loading  and  the  stiffness  and  volume  change  behavior  in  drained 
triaxial  compression.  Although  in  all  cases  studied  in  this  investigation 
higher  strengths  under  the  triaxial  stress  states  used  were  associated  with 
a greater  randomness  in  particle  long  axis  orientations,  the  relations 
between  strength  and  interparticle  contact  orientations  were  reversed  for 
samples  at  50  percent  and  80  percent  relative  density.  It  should  be  noted 
that  these  findings  are  somewhat  at  variance  with  those  for  the  lass 
extensive  and  preliminary  fabric-strength  study  reported  by  Mulilis  et  al 
(1975). 
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There  are  other  aspects  of  sample  preparation  methods  and  the 
fabric  differences  studied  here  that  could  contribute  to  the  observed 
differences  in  strength  and  volume  change  behavior.  Several  observations 
were  made  during  the  study  of  thin  sections  that  may  be  potentially 
significant.  Insufficient  study  of  them  has  been  made  to  allow  detailed 
quantitative  analysis.  They  are 

1.  Of  the  order  of  0.5  to  1.0  percent  of  the  grains  in  samples 
prepared  by  moist  tamping  and  moist  vibration  appeared  cracked 
or  fractured.  A schematic  drawing  of  the  nature  of  this 
observation  is  shown  in  Fig.  34.  A detailed  macroscopic  study 
of  grains  before  and  after  compaction  would  be  required  to 
confirm  that  cracking  resulted  from  compaction  stresses. 

2.  Larger  voids  were  observed  in  the  horizontal  thin  sections 
than  in  vertical  sections.  Figure  35  is  a schematic  diagram  of 
an  isolated  large  pore. 

3.  Grain  and  contact  orientations  are  less  easily  determined  in 
dense  them  in  loose  regions. 

There  are  other  fabric-related  factors  and  observations  that  have  not 
yet  been  studied  in  detail,  but  which  may  have  a major  influence  on  the 
liquefaction  behavior  of  sand.  They  are: 

1.  Variable  void  ratios  within  samples.  The  influences  of  both 
large  and  small  inhomogeneities  within  samples  in  relation  to  a 
given  average  relative  density  remain  largely  unknown. 

2.  Segregation  of  different  particle  sizes  within  the  sample. 

3.  Influences  of  stress  history,  including  overconsolidation,  and 


types  of  stress  system.  Triaxial  loading  was  used  for  this 
investigation.  Fabric-liquefaction  relationships  for  deformation 
in  simple  shear  may  not  be  the  same. 


4.  Different  particle  cluster  arrangements  as  a result  of  dif- 
ferent preparation  methods.  Pore  size  distribution  measure- 
ments might  provide  information  relative  to  this,  as  well  as 
relative  to  factors  1 and  2 above . 

5.  These  studies  as  well  as  those  by  Mulilis  et  al  (1975)  and 
Ladd  (1975)*  show  that  wet  preparation  methods  consistently 
give  higher  strengths  than  do  dry  preparation  methods.  Both 
moist  tamping  and  moist  vibration  give  higher  strength  than  do 
dry  tamping  or  dry  vibration  to  the  same  density.  Ladd's  studies 
show  that  the  rate  of  strain  development  with  increased  numbers 
of  loading  cycles  was  the  same  for  samples  of  a uniform  sand 
prepared  by  all  methods  to  a relative  density  of  60  percent.  On 
the  other  hand  the  rate  of  strain  development  was  quite  different 
for  dry  fabricated  and  moist  fabricated  samples  of  a more  well 
graded  sand  prepared  to  a relative  density  of  80  percent. 

6.  Chemical  effects  at  interparticle  contacts  that  may  develop  in 
moist  samples  owing  to  abrasion  during  tamping  or  vibration  and 
subsequent  consolidation  may  be  of  major  importance.  Chemical 
analysis  and  high  resolution  electron  microscopy  of  particle 
surfaces  before  and  after  compaction  might  provide  useful  data  on 
chemical  and  surface  texture  phenomena. 

7.  The  influence  of  sand  type  to  include  particle  shape  and 
gradation.  The  Monterey  No.  0 sand  used  for  this  study  is 
uniform  and  composed  of  rounded  grains  which  are  nearly  equi- 
dimensional. 

Finally,  it  is  clear  that  careful  study  is  needed  of  the  fabrics 
of  naturally  occurring  sand  deposits  and  sand  fills  in  the  field  if 

*Ladd,  R.  S.  (1975),  personal  communication. 
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proper  characterization  of  their  behavior  by  tests  on  reconstituted 
samples  in  the  laboratory  is  to  be  achieved.  In-situ  injection  of  sands 
i.i  the  field  to  stabilize  the  fabric  prior  to  sampling  should  enable 
recovery  of  samples  suitable  for  detailed  study. 


X.  CONCLUSION 


This  study  has  established  that  the  undrained  strength  of  Monterey 
No.  0 sand  under  cyclic  triaxial  loading  and  the  drained  triaxial  stress- 
strain  and  volume  change,  strain  behavior  bear  specific  relationships  to 
the  sand  fabric,  which  in  turn  is  determined  by  the  method  of  sample 
preparation.  Of  three  sample  preparation  methods  investigated  moist 
vibration  gave  the  most  random  orientations  of  particle  long  axes.  Dry 
pluviation  resulted  in  distinct  preferred  orientation  of  particle  long 
axes  in  the  horizontal  direction.  Moist  tamping  gave  intermediate  values. 
The  highest  strengths  were  associated  with  the  most  random  particle 
orientations. 

The  distribution  of  interparticle  contact  plane  orientations  appears 
to  be  an  important  factor  influencing  the  resistance  of  a sand  to  lique- 
faction. However,  the  influence  of  contact  orientation  distributions 
appears  to  vary  depending  on  relative  density. 

A clea1'  relationship  exists  between  the  stiffness  and  volume  change 
characteristics  in  drained  triaxial  compression  and  the  resistance  to 
liquefaction  under  undrained  triaxial  cyclic  loading.  This  correlation  is 
stronger  at  a relative  density  of  50  percent  than  at  a relative  density  of 
80  percent. 

Additional  factors  such  as  pore  size  and  particle  cluster  size  dis- 
tribution, particle  breakage  during  compaction,  particle  size  segregation 
within  samples,  and  chemical  effects  at  interparticld  contacts  need  further 
study  in  relation  to  their  influences  on  strength  and  liquefaction  behavior. 

Direct  application  of  the  results  of  this  study  to  the  selection  of 
the  appropriate  method  for  preparation  of  reconstituted  samples  must  await 
the  analysis  of  the  fabrics  of  undisturbed  samples  from  the  field. 
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FIG.  2 GRAIN  SIZE  DISTRIBUTION  OF  MONTEREY  NO.  0 SAND 
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(a)  Locations  of  Tnitio!  Cuts 


(b)  Horizontal  Section 


(c)  Vertico / Sections 

FIG.  4 THIN  SECTION  LOCATIONS  FROM  WITHIN 
RESIN-IMPREGNATED  SAMPLES 


FIG.  5 PHOTOMICROGRAPH  OF  VERTICAL  THIN  SECTION  FROM  SAMPLE  PREPARED  BY  PLUVIATION 
(50%  RELATIVE  DENSITY) 


FIG.  6 PHOTOMICROGRAPH  OF  VERTICAL  THIN  SECTION  EROM  SAMPLE  PREPARED  BE  MOIST  TAMPING 
(50%  RELATIVE  DENSITY) 


FIG.  7 PHOTOMICROGRAPH  OF  VERTICAL  THIN  SECTION  FROM  SAMPLE  PREPARED  BY  MOIST  VIBRATION 
(50%  RELATIVE  DENSITY) 


FIG.  8 PHOTOMICROGRAPH  OF  VERTICAL  THIN  SECTION  FROM  SAMPLE  PREPARED  BY  PLUVIATION 
(80%  RELATIVE  DENSITY) 


Pluviated  D 


Moist  Tamped  D 


FIG.  11  X-RADIOGRAPHS  OF  LONGITUDINAL  SECTIONS  FROM  SAMPLES  AT 
80%  RELATIVE  DENSITY 
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FIG.  12  UNIVERSAL  STAGE  LOOKING  DOWN  THE 
VERTICAL  ROTATION  AXIS 


FIG.  13 
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MEASUREMENT  OF  NORMALS  (Nj.  , Ni  ) 
TO  TANGENT  PLANE  ( tt  . ) 


p/uviotion 


FIG.  15  LIQUEFACTION  BEHAVIOR  OF  MONTEREY  NO.  0 SAND  PREPARED  TO  A RELATIVE  DENSITY 
OF  80%  BY  PLUVIATION  AND  MOIST  TAMPING 
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FIG.  18  INFLUENCE  OF  COMPACTION  METHOD  ON  THE  DRAINED 
TRIAXIAL  COMPRESSION  BEHAVIOR  OF  MONTEREY  NO. 
SAND  AT  80%  RELATIVE  DENSITY 


INFLUENCE  OF  PLUVIATION  NOZZLE  SIZE  ON  THE  DRAINED 
COMPRESSION  BEHAVIOR  OF  MONTEREY  NO.  0 SAND  AT  80% 
RELATIVE  DENSITY 
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FIG.  20  HISTOGRAMS  OF  PARTICLE  LONG  AXIS  ORIENTATIONS  FOR  SAMPLES  OF 

MONTEREY  NO.  0 SAND  PREPARED  TO  50%  RELATIVE  DENSITY  BY  DIFFERENT 


(o)  Pluvioted  L =18  3%,  6*86 


ROSE  DIAGRAMS  OF  PARTICLE  LONG  AXIS 
ORIENTATIONS  FOR  SAMPLES  OF  MONTERE 
SAND  PREPARED  TO  50%  RELATIVE  DENSI 
DIFFERENT  METHODS 
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FIG.  23  ROSE  DIAGRAMS  OF  PARTICLE  LONG  AXIS 

ORIENTATIONS  FOR  SAMPLES  OF  MONTEREY  NO.  0 
SAND  PREPARED  TO  80%  RELATIVE  DENSITY  BY 
DIFFERENT  METHODS 


FIG.  28  EQUAL  AREA  STEREONET  SHOWING  DISTRIBUTION  OF 
INTERPARTICLE  CONTACT  NORMALS  IN  A SAMPLE 
PREPARED  BY  MOIST  TAMPING  TO  A RELATIVE  DENSITY 
OF  80% 


FIG.  29  INTERPARTICLE  CONTACT  NORMAL  DISTRIBUTION  FUNCTIONS  FOR  SAMPLES  PREPARED 
BY  PLUVIAL  COMPACTION 
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INTERPARTICLE  CONTACT  NORMAL  DISTRIBUTIONS  FOR  SAMPLES  PREPARED  TO 
RELATIVE  DENSITY  OF  50%  BY  THREE  METHODS 
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FIG.  33  INTERPARTICLE  CONTACT  NORMAL  DISTRIBUTIONS  FOR  SAMPLES  PREPARED  TO  A 
RELATIVE  DENSITY  OF  80%  BY  TWO  METHODS 


FIG.  34  SCHEMATIC  DIAGRAM  OF  TENSION 
CRACKS  IN  GRAINS  WITHIN  MOIST 
VIBRATED  AND  MOIST  TAMPED 
SAMPLES  OF  MONTEREY  NO.  0 SAND 


FIG.  35 


SCHEMATIC  DIAGRAM  OF  A LARGE 
PORE  AS  SEEN  IN  HORIZONTAL 
THIN  SECTION  (~60X  MAGNIFICATION) 
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In  accordance  with  TR  70-2-3,  paragraph  6c(l)(b), 
dated  15  February  1973,  a facsimile  catalog  card 
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